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ABSTRACT 

The aim of this work is to identify the physical processes that occur in the 
network and contribute to its dynamics and heating. We model the network as 
consisting of individual flux tubes with a non-potential field structure that are 
located in intergranular lanes. With a typical horizontal size of 200 km at the 
base of the photosphere, they expand upward and merge with their neighbors 
at a height of about 600 km. Above a height of approximately 1000 km the 
magnetic held starts to become uniform. 

Waves are generated in this medium by means of motions at the lower bound¬ 
ary. We focus on transverse driving, which generates both fast and slow waves 
within a hux tube and acoustic waves at the interface of the tube and the held- 
free medium. The acoustic waves at the interface are due to compression of the 
gas on one side of the hux tube and expansion on the other. These waves travel 
upward along the two sides of the (2D) hux tube and enter it, where they be¬ 
come longitudinal waves. For impulsive excitation with a time constant of 120 s, 
we hnd that a dominant feature is the creation of vortical motions that prop¬ 
agate upward. We have identihed an efficient mechanism for the generation of 
longitudinal waves and shock formation in the chromospheric part of the hux 
concentration. We examine some broad implications of our results. 

Subject headings: MHD — Sun: chromosphere — Sun: magnetic helds —- Sun: 
oscillations 
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1. Introduction 

The quiet solar chromosphere is bifurcated into the magnetic network on the boundary 
of supergranulation cells and the largely held-free internetwork medium in the cell interior, 
with respective hlling factors in the ratio of 2 to 3 in the middle chromosphere (i.e., at a 
height of about 1 Mm above the level of unit optical depth at 5,000 A). The network is about 
30% brighter in Ca II K 2 v emission than the internet chromosphere (Skumanich, Smythe & 
Frazier 1975). Both the network and the internet medium show bright points (BPs), which 
are prominent in the emission peaks in the cores of the Ca II H and K lines, formed in 
the middle chromosphere. However, the dynamical and spectral properties of network and 
internet BPs are quite different. In internetwork areas the chromospheric velocity power 
spectrum is dominated by oscillations with frequencies at and above the acoustic cutoff 
frequency (period 3 min in upper photosphere, e.g., Rutten & Uitenbroek 1991), whereas 
in the network, Ca II H line center velocity and intensity power spectra are dominated by 
low-frequency oscillations with periods of 7-20 min (Bites, Rutten & Kalkofen 1993). These 
long-period waves have also been observed at larger heights (Curdt & Heinzel 1998; McAteer 
et al. 2002, 2003). Furthermore, network BPs show high emission in the line core most of 
the time (Bites et al. 1993), whereas internetnetwork locations have a bright phase only 
about 5% to 10 % of the time (von Uexkiill & Kneer 1995). Finally, 80% of network BPs 
have symmetric prohles in the line core, with more or less equal intensity in the blue and 
red emission peaks on either side of the central absorption. In contrast, only about 30% of 
internet BPs have symmetric prohles, and 40% have blue-peak enhancements (Grossmann- 
Doerth, Kneer & v.Uexkiill 1974); hence the name H 2 v or K 2 v bright point. 

While the qualitative properties of internetwork BPs are reasonably well understood, 
including their formation in upward-propagating acoustic shocks that encounter downward- 
howing gas (Carlsson & Stein 1995), this is not true for network BPs. The physical processes 
that heat the magnetic network have not been fully identihed. Are network BPs heated by 
wave dissipation, and if so, what is the nature of these waves? How can we understand 
the observed relatively constant emission and symmetric line prohles of network BPs? The 
source of energy for network BPs is likely to be magnetohydrodynamic waves. Possible 
candidates are: 

• Kink (transverse) waves, generated inside hux tubes by the buheting action of granules; 

• Bongitudinal waves, generated by pressure huctuations inside hux tubes; 

• Torsional (Alfven) waves, generated inside hux tubes; 

• acoustic waves, generated in the held-free atmosphere surrounding hux tubes, that 



penetrate into the tubes; 

• Acoustic waves, generated at the interface of flux tubes and the outside medium, that 
also penetrate into the flux tubes. 

The first three are well-known flux tube modes (e.g.. Spruit 1981). Several investigations 
have focused on the generation and propagation of transverse and longitudinal wave modes 
and their dissipation in the chromosphere (e.g., Zhugzhda, Bromm & Ulmschneider 1995; 
Fawzy et ah 2002; Ulmschneider 2003 and references therein). Torsional waves have received 
some attention (e.g., Hollweg 1982; Noble, Musielak & Ulmschneider 2003). The fourth 
wave type is an obvious source since it is responsible for heating the internetwork medium. 
However, the importance of the hfth wave type, which we found during the course of this 
investigation, does not appear to have been adequately recognized earlier. It will be discussed 
in more detail later on. 

The present study is a continuation of earlier work on the excitation of transverse and 
longitudinal waves in magnetic flux tubes by the impact of fast granules on flux tubes (Hasan 
& Kalkofen 1999), as observed by Muller & Roudier (1992) and Muller et ah (1994), and 
following the investigation by Choudhuri, Auffret & Priest (1993), who studied the generation 
of kink waves by footpoint motion of flux tubes. The observational signature of the modeled 
process was highly intermittent in radiation emerging in the H and K lines, contrary to 
observations. By adding waves that were generated by high-frequency motion due to the 
turbulence of the medium surrounding flux tubes (Hasan, Kalkofen & van Ballegooijen 2000), 
the energy injection into the gas inside a flux tube became less intermittent, and the time 
variation of the emergent radiation was in better agreement with the more steady observed 
intensity from the magnetic network. 

The above studies modeled wave excitation and propagation in terms of the Klein- 
Gordon equation, motivated by the identihcation of the power peak near 7 min in the ob¬ 
served power spectrum (Bites et al. 1993; Curdt & Heinzel 1998; McAteer et ah 2002) with 
the cutoff period of kink waves in thin magnetic flux tubes (Kalkofen 1997). This analysis 
was based on a linear approximation in which the longitudinal and transverse waves are 
decoupled. However, it is well known that the velocity amplitude v{z) for the two modes 
increases with height (for an isothermal atmosphere, as v oc exp[z/4H], where H is the 
pressure scale height), so that the motions are expected to become supersonic higher up in 
the atmosphere. At such heights, nonlinear effects become important, leading to a coupling 
between the transverse and longitudinal modes. Some progress on this question has been 
made in one dimension, using the nonlinear equations for a thin flux tube, by Ulmschnei¬ 
der, Zahringer & Musielak (1991), Huang, Musielak & Ulmschneider (1995), and Zhugzhda, 



Bromm & Ulmschneider (1995), and more recently by Hasan et al. (2003) and Hasan & Ulm- 
schneider (2004), who examined mode conpling between transverse and longitndinal modes 
in the magnetic network. By solving the nonlinear, time-dependent MHD eqnations it was 
fonnd that signihcant longitndinal wave generation occurs in the photosphere, typically for 
Mach numbers as low as 0.2, and that the onset of shock formation occurs at heights of 
about 600 km above the photospheric base, accompanied by heating (Hasan et al. 2003, 
Hasan & Ulmschneider 2004). The efficiency of mode coupling was found to depend on the 
magnetic held strength in the network and is a maximum for held strengths corresponding to 
(3 ^ 0.2, when the kink and tube wave speeds are almost identical. This can have interesting 
observational implications. Furthermore, even when the two speeds are diherent, once shock 
formation occurs, the longitudinal and transverse shocks exhibit strong mode coupling. 

The above studies on the magnetic network make use of two important idealizations: 
they assume that the magnetic hux tubes are thin, an approximation that becomes invalid 
at about the height of formation of the emission peaks in the cores of the H and K lines; and 
they neglect the interaction of neighboring flux tubes. Some progress in this direction has 
recently been made by Rosenthal et al. (2002) and Bogdan et al. (2003), who studied wave 
propagation in a two-dimensional stratihed atmosphere, assuming a potential magnetic held 
to model the network and internetwork regions on the Sun. They examined the propagation 
of waves that are excited from a spatially localized source in the photosphere. Their results 
indicate that there is strong mode coupling between fast and slow waves at the so-called 
magnetic canopy, which they identify with regions where the magnetic and gas pressures are 
comparable. As a consequence of the potential-held approximation, some magnetic held lines 
are nearly horizontal even at the base of the held. Such a model may not be appropriate for 
a network patch, which is perhaps better idealized as a collection of vertical tubes (Cranmer 
& van Ballegooijen 2005). 

Thus the problem that we address concerns wave propagation in regions that are largely 
representative of individual structures in the magnetic network — this is diherent from the 
one analysed by Rosenthal et al. and Bogdan et al. Our initial conhguration consists of 
hux tubes in 2-D magnetostatic equilibrium with a sharp interface between the tube and 
the surrounding gas. Waves are generated in this medium by means of transverse motion 
at the lower boundary, which displaces the entire hux tube, unlike the problem studied in 
the above papers where the source region is conhned to a portion of the magnetic structure. 
This can have interesting consequences, some of which were unrecognized thus far. These 
will be explored in later sections of the paper. Our calculations indicate the presence of a 
new and efficient mechanism for longitudinal wave generation and shock formation in the 
chromosphere. 
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The present study forms the hrst of a series devoted to a detailed investigation of wave 
propagation in the magnetic network. As a hrst step, we go beyond the thin hux tube ap¬ 
proximation and employ a two-dimensional treatment in slab geometry, similar to Rosenthal 
et ah and Bogdan et al. It is not possible (within the 2-D framework) to examine kink and 
sausage waves, which would require a three-dimensional treatment, but, nevertheless, there 
are still many important physical effects that this approach allows us to examine. For sim¬ 
plicity, we neglect non-adiabatic effects at this stage, which we will include in a separate 
paper. 

Several papers have looked at MHD waves in various geometries in the solar atmosphere 
using multi-dimensional simulations (e.g., Shibata 1983; Cargill, Spicer, & Zalesak 1997; 
Ofman & Davila 1998; Sakai, Igrarashi & Kawata 1998; Ofman, Nakariakov, & DeForest 
1999; additional references can be found in Bogdan et al. 2003). 

The organization of this paper is as follows: in Sect. 2 we describe the initial equilibrium 
model, its construction and in Sect. 3 the numerical method of solution along with the driving 
mechanism. The results of our calculation are presented in Sect. 4 followed by a discussion 
and summary in Sect. 5. 


2. Model 

Following Cranmer & van Ballegooijen (2005), we treat a network element (typical flux 
~ 3 X 10^® Mx) to consist of a collection of smaller flux tubes that are spatially separated 
from one another in the photosphere. The gas pressure in the atmosphere decreases with 
increasing height, causing a lateral expansion of the flux tubes. Neighboring flux tubes within 
the network element merge into a monolithic structure at some height. Above this “merging 
height” the network element consists of a single thick flux tube that further expands with 
height. The outer edge of this tube forms a magnetic canopy that overlies the neighboring 
supergranular cells. A second merging occurs when neighboring network elements come 
together at a canopy height. Figure 1 schematically shows the picture we have for the 
network held structure. It consists of three distinct regions: 

1. Photospheric region up to about 0.6 Mm, consisting of individual hux tubes, with a 
typical diameter of 200 km in the low photosphere. Their foot points are located in 
intergranular lanes and separated from one another by about the diameter of a granule 
(~ 1 Mm). They expand upward and merge with their neighbors at a height of about 
0.6 Mm. 

2. Lower chromosphere, between the heights of 0.6 and 1 Mm, the merged network hux 
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element expands laterally over the surrounding supergranular cell-center and overlying 
field-free chromosphere; 

3. Upper chromosphere and corona, between 1 and 12 Mm, the fully merged magnetic field 
fills the available volume. At larger heights, the field expands primarily in the vertical 
direction and becomes more or less uniform. However, at lower heights, between 1 
and 2 Mm, the field strength varies significantly with horizontal position, and the field 
strength directly above the flux tubes (left side of Figure 1) is much larger than above 
the supergranular cell center (right side of Figure 1). 

Our model is based on the idea that the base of flux tubes is located in subsurface layers 
where convective flows may be different from those in the visible photosphere. Flux tubes 
occur in regions with convective downflows below intergranular lanes. These downflows are 
likely to be highly turbulent, involving lateral motions that produce transverse waves in 
the flux tubes. When the upward propagating waves reach the photosphere, they cause 
horizontal motions of flux tubes relative to their local surroundings. This generates excess 
pressure on the leading edge of a flux tube, and a pressure deficit on the trailing edge. Our 
two-dimensional MHD calculations (see section Sect. 3) indicate that these pressure pulses 
produce an upward-velocity pulse on the leading side of the flux tube and a downward pulse 
on the trailing side. These pressure pulses and vertical flows are an integral part of the MHD 
wave, so the transverse and longitudinal motions are strongly coupled. 


2.1. Initial Two-dimensional Magnetostatic Model 

Let us consider an individual flux tube at the base of a magnetic network element on the 
quiet sun; the region of interest is indicated by the small box in Figure 1. At heights below 
about 600 km, the flux tubes are spatially distinct from one another, and are embedded in a 
field-free “external” medium. At ~ 1000 km, the flux tubes merge into a more uniform field 
(upper part of small box). The structure of the flux tube at the initial instant is assumed to 
be in static equilibrium and is determined by the magnetostatic force balance equation: 

-Vp -F pg -F ^(V X B) X B = 0, (1) 

Itt 

where g = —gz is the gravitational acceleration, p is the gas pressure, p is the density, and 
B is the magnetic field. The third term in (1) describes the Lorentz force due to electric 
currents at the boundary between the flux tube and its local surroundings. 

We have developed a numerical code for solving equation (1) in two dimensions; all 



7 


quantities are assumed to be independent of the horizontal y coordinate, so the flux tubes 
are approximated by sheets. We consider a rectangular domain x = [0, L] and 2 : = [0,77], 
representing one half of a flux sheet; a; = 0 is the flux sheet axis, x = L is the interface with 
the neighboring sheet on the right, and = 0 is the base of the photosphere. The magnetic 
field is written in terms of a flux function A{x, z): 


Bx 


dz' ^ dx 


( 2 ) 


The boundary conditions are ^4(0, z) = 0 and A{L, z) = Amax, where is one half of the 
total magnetic flux within the sheet. The gas pressure and density are written as: 


P = Pm,(2) [1+(3„'F(/l)] , (3) 

P = Pi„,W[l+^(„-‘F(^)], (4) 


where Pint(^) is the internal gas pressure as function of height along the axis of the flux sheet; 
Pint( 2 :) is the internal density, Pint{z) = —g~^dpira./dz] Pq is the ratio of gas- and magnetic 
pressures at the base (on axis); and F{A) is a function describing the variation of gas pressure 
across held lines. Note that the temperature depends only on height {p/p is independent 
of x). The function F{A) varies from zero on the axis of the sheet to F(Aniax) = 1 in the 
external medium, resulting in distributed electric currents at the interface between the sheet 
and its local surroundings. 


Inserting (2), (3) and (4) into (1) yields 


VM 4npi^t{z)Po ^ 


dF 

dA 


0 , 


which can be solved by minimizing the following Lagrangian: 


W 



47rpint(^)/7o ^F(A) 


dx dz. 


( 5 ) 

( 6 ) 


The minimization is done by varying the A-values on a grid of 120 by 240 cells, using 
the conjugate-gradient method (Press et al. 1992). A similar technique was used in (Hasan 
et al. 2003), where we constructed a model of a very thin hux tube (see the Appendix of this 
paper), and in Cranmer & van Ballegooijen (2005). In these papers a Lagrangian description 
was used, and the radial positions r{^,z) varied on a fixed grid of <h (the flux function in 
cylindrical coordinates) and height ; 2 . In the present paper we vary A{x, z) on a hxed grid 
of X and 2 :. 

A model containing one whole flux sheet can be obtained by mirroring the held with 
respect to a; = 0, and models containing mnltiple sheets can be obtained by fnrther mirroring 



with respect to x = L. The resulting helds will be used as initial conditions for the two- 
dimensional MHD calculation (see next subsection). 


We perform calculations for a single sheet using the following form for F{A): 


F{A) 


I {A/A^^^f - \ {A/A^^^f if A < 

1 otherwise. 


(7) 


This produced a smooth transition of gas pressure from the interior to the exterior of the 
sheet, and the residual magnetic held in the external medium was very small (~ 1 Gauss). 
The internal pressure as function of height was approximated as a sum of two exponentials: 


Pint(^) =Piexp(-^/hfi) -Fp2exp(-2:/if2), 


( 8 ) 


with a photospheric pressure scale height of Fli = 110 km and a chromospheric scale height 
of H 2 = 220 km. To obtain kilogauss helds in the photosphere, we used Pq = 0.5, so that 
the external gas pressure is three times the internal gas pressure. 

In the present calculations we consider a full hux sheet which is placed symmetrically 
with respect to the midpoint of the computational domain, which we take to be a square 
of size 1200 km with a uniform grid of 240 x 240 cells, corresponding to a mesh spacing of 
5 km in either direction. The axis of the hux tube is now located at a; = 600 km. The width 
of the sheet at = 0 is approximately 200 km, and the central held strength is 1530 G. 

Table 1 provides some of the parameters of the equilibrium model on the tube axis and 
in the ambient medium at the base {z = 0 ) and the top {z = 1200 km) of the computational 
domain. Figures 2(a) and (b) depict the variation with height of the temperature and 
pressure, respectively. It should be noted that the temperature by assumption is constant 
in the horizontal direction. The sound speed C 5 corresponding to the temperature shown in 
Figure 2(a) varies from 7.2 km s“^ at = 0 to 9.3 km s“^ aX z = 1200 km. In Figure 2(b), 
the pressure is shown on diherent held lines corresponding to a fractional hux / = AjA^aax 
as measured from the left edge of the computational domain. Thus, / = 0.5 corresponds to 
the tube axis. As expected, the gas pressure on each held line increases as one approaches 
the tube boundary. The dashed curve (/ = 0.95) essentially depicts the height variation of 
the pressure in the ambient medium. 

Figures 3(a) and (b) show the height variation of the magnetic held strength B and of 
P on various held lines, parameterized, by / = 0.5, 0.8, 0.9. The held strength on each of the 
held lines drops rapidly with in the hrst few hundred kilometers after which it approaches 
a constant value of about 100 G. On the other hand, P is practically constant with in the 
lower region of the atmosphere, which is similar to the behavior one hnds in an thin hux tube 
with equal internal and external temperature at the same height (e.g. see Hasan et ah 2003), 
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and where both and p have the same height dependence. In the npper part of the flnx 
tnbe, B is constant and (3 drops off sharply with essentially mimicking the p dependence. 
The Alfven speed, which is related to f3 by the approximate relation va = cs/\/]3^ is almost 
constant (with height) in the lower regions of the atmosphere in the tnbe, bnt increases 
sharply with in the higher layers. 

Fignre 4 shows the variation with horizontal distance x of the vertical component of 
the magnetic held Bz at varions heights =0 (black solid line), 500 km (black dotted line) 
and 1000 km (black dashed line). The red cnrves show the corresponding variation of B^ 
measnred with respect to the scale on the right. The solid lines (black and red) clearly 
show that the magnetic held at the base is conhned to a width of 200 km. At a height 
of 1000 km, Bx becomes very small (a few Ganss) and the held is essentially vertical and 
constant (approximately 100 G) in the horizontal direction. 


3. Method and Boundary Conditions 

We consider wave generation in the conhgnration described in the previons section 
by perturbing the lower boundary and solving the two-dimensional magnetohydrodynamic 
(MHD) equations in conservation form for an inviscid adiabatic huid. These consist of the 
usual continuity, momentum, entropy (without sources) and magnetic induction equations 
(for details see Steiner, Knolker & Schiissler 1994). The unknown variables are the density, 
momentum, entropy per unit mass and the magnetic held. We assume that the plasma 
consists of fully ionized hydrogen with a mean molecular weight of 1.297. The temperature 
is computed from the specihc entropy and the pressure is found using the ideal gas law. 

The above MHD equations are solved following the numerical procedure given by Steiner 
et al. (1994). Briehy, the equations are discretized on a two-dimensional mesh using a hnite- 
volume method which has the advantage of preserving div B = 0 to machine accuracy. The 
method employs hnite diherences to compute numerical huxes based on the hux-corrected 
transport (FGT) scheme of Oran and Boris (1987). The time integration is explicit and has 
second-order accuracy in the time step. Small time steps are required to satisfy the Gourant 
condition in the upper part of the domain, where the Alfven speed is large. 

Periodic boundary conditions are used at the horizontal boundaries. At the top bound¬ 
ary, (a) the vertical and horizontal components of momentum are set to zero; (b) the density 
is determined using linear logarithmic extrapolation; (c) the horizontal component of the 
magnetic held and temperature are set equal to the corresponding values at the preceding 
interior point. The vertical component of the magnetic held is determined using the con- 
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dition div il=0. Similar conditions are used at the lower boundary, except for the density, 
temperature and horizontal component of the velocity (or momentum). The density is ob¬ 
tained using cubic spline extrapolation, the temperature is kept constant at its initial value 
and Vx sX z = Q specihed as follows: 

na;(a;, 0,f) = nosin(27rt/P), (9) 

where vq denotes the amplitude of the horizontal motion and P the wave period. This 
form was chosen to simulate the effect of transverse motion of the lower boundary. For 
simplicity, we assume that all points at the lower boundary have this motion, since this does 
not generate any waves in the ambient medium, other than at the interface with the flux 
tube, as we shall see later on. 


4. Dynamics of a flux sheet 

In Sect. 2, we presented a model for a single flux sheet in static equilibrium. The 
stability of this equilibrium was checked by solving the time-dependent two-dimensional 
MHD equations without any external driving (assuming rigid boundaries in the vertical 
direction). We found that this equilibrium was maintained to high accuracy. The maximum 
amplitude of the flows was no more than a few meters per second over time scales greater 
than the sound travel time (in both the vertical and horizontal directions), corresponding to 
over 10000 time steps. 

Let us now consider wave generation in the equilibrium atmosphere by means of a 
transverse motion at the lower boundary (z = 0) of the atmosphere, which displaces the flux 
sheet. We focus on two limiting cases, corresponding to impulsive and periodic excitation, 
respectively. 


4.1. Impulsive excitatiou 

This case corresponds to uniform displacement of the lower boundary at = 0 to the 
right which lasts half a wave period (P = 240 s) and then stops. The peak transverse 
velocity is vq = 750 m s“L Figures 5(a)-(d) show, at times of 52 s, 82 s, 109 s and 136 s, 
respectively, the velocity field (arrows), the magnetic field (black lines) and the temperature 
change AT = T — Tq, where Tq is the initial temperature at each height. The maximum 
value of the velocity and the color table for the temperature are shown on the right of each 
hgure. We have omitted velocities smaller than 30 m s“^. The white lines denote contours 
of constant beta for values of (3 = 0.1, 0.5, 1.0 (thick line) and 10. 



The horizontal motion of the flux tube at the lower boundary pushes the held lines 
uniformly to the right. In the hux tube, this motion is communicated to the upper layers as 
a fast mode, which travels with the local Alfven speed. The transverse motion of the held 
lines compresses the gas on the right interface of the tube and the ambient medium (above 
the base, where there is no horizontal motion in the held free medium). 

Figure 6 clearly shows at t = 82 s the development of a pressure enhancement and 
depletion, respectively, at the right and left interfaces of the tube with respect to the held 
free gas. This essentially creates a pressure dipole which in turn generates a vortex motion 
with uphow and downhow motions on the right and left sides, respectively. The top of the 
vortex motion is in the opposite direction to that at the base. From Figure 5 we hud that, 
as time proceeds, the vortex grows in size and also moves upward. This motion carries the 
held lines, pushing them slightly counterclockwise in the lower regions of atmosphere. 

As the vortex extends upwards, there is a hint of a shock forming at its upper edge, 
as can be seen in Figure 5(c). Once the transverse driving motion stops at the base, the 
vortex motions diminish and the how is now guided along the held lines. This is evident 
in Figure 5(d) which shows at t = 136 s negligible hows near the base of the hux tube 
and large hows almost aligned with the magnetic held in the upper atmosphere. The large 
uphows generate a shock which can be discerned at a height of about 900 km along with a 
temperature enhancement greater than 400 K. 

Figures 7(a) and (b) denote the 2 : variation of Vg and Vn, the held-aligned and normal 
velocity components respectively at the times of 82, 109 and 150 s along a specihc held 
line. We choose a held line parameterized by / = 0.6, which is just to the right of the 
hux tube axis. Let us hrst consider the behavior of the held-aligned how, shown in Figure 
7(a). As time proceeds, the vertical component of the velocity pulse increases with height in 
magnitude on account of the density stratihcation and it steepens on account of nonlinear 
ehects as it propagates upward. At t = 150 s the pulse resembles a shock. The presence of a 
wake behind the shock should also be noted, which is a generic feature associated with the 
propagation of a pulse in a vertically stratihed atmosphere (e.g. Hasan et al. 2003). Figure 
7(b) shows the variation of the normal component of the how. The curves corresponding 
to t = 82 and 109 s show the normal component prohles during the impulsive phase (i.e., 
during the interval of transverse driving motion at the base is present). The dashed curve 
shows the prohle at t = 150 s, after the transverse motion at the base has halted. This prohle 
is broadened due to the fact that the Alfven speed, which is the characteristic propagation 
speed of the transverse pulse, increases sharply with height. In the lower regions of the 
atmosphere, where (3 ~ 0.6, va ~ 10 km s“^, whereas at a height of 600 km, it has increased 
to about 26 km s“^. 



4.2. Periodic excitation 


We now consider the periodic excitation of waves due to transverse driving of the lower 
boundary, with a period of 24 s (similar to Rosenthal et ah and Bogdan et ah), though with 
a higher amplitude of 750 m s“^. Figures 8(a)-(d) show the wave pattern that develops in 
this case. As before, the colors are used to denote the temperature perturbation with respect 
to the initial value at each height and the white curves depict contours of constant 13. 

The horizontal motion of the tube at the lower boundary is a source of acoustic waves just 
outside the flux tube boundary in the ambient medium, in which they propagate isotropically 
and near the interface as ‘spherical’ acoustic waves. The yellow and blue strips associated 
with AT, with an almost constant separation, clearly show the acoustic waves propagating 
outwards at the sound speed from their source near the bottom edge of the tube. The 
wavelength of these waves is approximately constant since the acoustic speed varies weakly 
in the lower part of the computational domain. 

In the flux tube, the horizontal driving motions generate both fast and slow modes. 
Close to the tube axis, the held is strong (/5 < 1) and the transverse motion generates a fast 
mode which propagates upwards at the Alfven speed. Its propagation can be clearly seen in 
Figure 9, which depicts the variation of the normal velocity component Vn as a color contour 
plot. The fast wave front shows an asymmetry in propagation since the Alfven speed is 
largest on the tube axis and decreases in the horizontal direction. The increase in separation 
in the vertical direction between the color peaks clearly shows the increase in wavelength 
due to the increase in Alfven speed with ;2, which increases very rapidly with height from 
about 10 km s“^ at the base to about 100 km s“^ at the top of our computational box. 

The horizontal motions in the central part of the flux tube produce compressions and 
decompressions of the gas at the edge of the tube, where the magnetic held is weaker (/5 > 1). 
The resulting pressure variations are periodic in height and time, and are 180° out of phase 
on the two sides of the hux tube. The associated pressure gradients drive periodic vertical 
hows that propagate upward along the edge of the hux tube at the speed of the transverse 
wave. These vertical hows are an integral part of the hux tube wave, so the wave has both 
transverse and longitudinal character in diherent parts of the tube. As the wave travels 
upward, it reaches a height where the magnetic held at the edge of the tube becomes more 
horizontal, and longitudinal motions are driven into the central part of the hux tube where 
/5 < 1; this process is analogous to the mode conversion discussed by Bogdan et al. (2003). 
In the low-/5 region, the wave consists of two almost decoupled parts: a fast mode that 
propagates in all directions, described in the preceding paragraph (see Figure 9), and two 
slow modes that are guided along the held lines. The slow modes on opposite sides of tube 
axis are 180° out of phase with each other, and their amplitudes increase with height due to 
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the stratification. Being mainly acoustic in nature, they produce a compression and heating 
of the plasma in the chromosphere. This can be discerned for instance in Figure 8(d), where 
the wave front located at a height of about a 1000 km produces heating greater than 500 
K. Figure 10 shows the variation of the held-aligned velocity component, Vg, on a held line 
/ = 0.6 at three different times, shown by the different curves. One can see the amplitude 
increase and steepening of the wave due to nonlinear effects as it propagates upward at the 
acoustic speed of around 8 km s“^. This clearly establishes the nature of this wave to be 
acoustic. It is interesting to note the occurrence of weak shocks at a Mach number of about 
0 . 2 . 


5. Discussion and Summary 

The present investigation is a continuation of our earlier work on the dynamics of the 
magnetic network on the Sun. In our previous studies we used a simplihed picture of mag¬ 
netic elements in terms of a quasi one-dimensional treatment based on the thin hux tube 
approximation. This approach was useful in providing a qualitative picture of wave propa¬ 
gation and nonlinear mode coupling in the lower regions of the atmosphere where the thin 
flux tube approximation is applicable. However, as already pointed out in Sect. 1, this ap¬ 
proximation breaks down in the chromosphere due to the expansion of the flux tube with 
height as well as due to the merging of different tubes. We overcome these limitations by 
using a two-dimensional treatment. In the first part of this study we focus on a single flux 
tube and examine the consequences when the lower boundary is perturbed by a transverse 
motion. Future work will extend the present analysis to multiple tubes as well as other types 
of excitation mechanisms. 

The choice of using a transverse velocity perturbation at the lower boundary was to some 
extent motivated by observations, particularly those of Muller and Rodier (1992) and Muller 
et al. (1994) who studied the footpoint motions of a large number of network bright points, 
generally regarded as a proxy for magnetic elements. We considered the driving motions to 
occur at a hxed height in the atmosphere. In reality, we expect that the displacement of the 
flux tubes occurs in response to turbulent motions below the photospheric base. These flows 
are absent or weak on the surface. Thus, the pressure fluctuations in our model are a direct 
consequence of the relative motion of the flux tube with respect to its local surroundings. 
These play a key role in driving the vortex motions that propagate high up in the flux tube. 
It appears reasonable to expect that such pressure fluctuations indeed occur in the solar 
network. 

We chose two limiting forms for the time behavior of the perturbation as an idealization 



-14- 


for drivers that (a) generate discrete pulses well separated in time and, (b) are periodic and 
continuous. The values of the time constant for the pulse in the first case and the wave 
period in the second case were chosen because of practical considerations related to the limit 
of not letting the simulation exceed a total time where reflections from the top boundary 
become important. 

We hnd that the transverse driving motions at the lower boundary lead to strong pres¬ 
sure perturbations in the field-free medium at the tube interface. This is an efficient process 
for generating vortical motions, recognized earlier by Shibata (1983) in two-dimensional sim¬ 
ulations, who found that such motions can arise due to a pressure perturbation applied at 
the base of a uniformly magnetized stratihed atmosphere (/? > 1). In our simulations the 
strong pressure fluctuations are essentially localized in regions where (3 is greater than unity. 
The observational signature associated with vortex formation would be the prediction of 
simultaneous upward and downward motions on opposite sides of flux tubes. Under what 
conditions does a vortex form? It appears that for wave periods that are sufficiently long 
(compared to the Alfven travel time) the dynamics is likely to be dominated by vortical 
motions, but it is difficult at this stage to say anything more dehnitive. 

An important consequence of our calculations is that the interface between the flux tube 
and the ambient medium is both a source of acoustic waves in the ambient medium as well 
as fast and slow waves just inside the flux tube. This can have interesting consequences: the 
acoustic waves that travel isotropically in the field free medium will impinge on neighboring 
flux tubes in the network and excite waves in them. On the other hand the modes generated 
in the high (3 region of the tube (at the interface) undergo mode conversion as they propagate 
upward and enter the region where ~ 1, as also found by Rosenthal et ah (2002) and 
Bogdan et ah (2003). In the upper part of the flux tube where the held is essentially vertical 
and /? 1, the acoustic-like longitudinal motions steepen and form shocks accompanied by 

heating, as also found out for instance by Hasan & Ulmschneider (2004) using a thin hux 
tube calculation. There is, however, an important qualitative difference between the present 
results and those of Zhugzhda et al. 1995, Hasan et al. (2003) and Hasan & Ulmschneider 
(2004). In the latter, longitudinal motions are generated from transverse driving motions at 
the hux tube footpoints as a consequence of nonlinear mode coupling that is most efficient 
when the transverse and longitudinal mode speeds are comparable (roughly for (3 ~ 0.2). 
In the present work, the mode conversion occurs essentially in the linear regime and is 
a consequence of the spatial variation of (3 in the tube. We hnd that transverse driving 
motions with velocities less than 1 km s“^ are sufficient to produce strong shocks in the 
chromosphere. 

The present work is in the spirit of the investigations by Rosenthal et al. (2002) and 



Bogdan et al. (2003) on wave propagation in two-dimensional magnetic structures. There are, 
however, qualitative differences between the initial state and the form of the perturbation 
used for driving the atmosphere in the two sets of calculations. For instance, the initial 
equilibrium conhguration adopted by the above authors assumes a potential held that is 
typical of the large-scale pattern connecting different network regions. Bogdan et al. (2003) 
consider a unipolar magnetic structure which at the base is 2000 km wide surrounded by 
smaller magnetized regions (about 750 km wide) of opposite polarity (as shown in Figure 
1 of their paper). At larger horizontal distance from the hux concentration, the held is 
essentially horizontal though much weaker. Thus, in the above conhguration the entire 
atmosphere is magnetized. On the other hand, we consider a discrete unipolar vertical 
hux, with a transverse dimension of about 200 km, in magnetostatic equilibrium (and non¬ 
potential held). Our equilibrium model, possessing a sharp interface between the tube and 
the ambient held-free medium, is perhaps more representative of a hux element in the solar 
network. A minor diherence is that the atmosphere in our model is not isothermal in the 
vertical direction but has a temperature increase with height similar to that in the solar 
chromosphere. 

Secondly, in our model, the whole hux tube at the photosphere is displaced and not 
just a small region inside it. The driving motion in the simulations of Bogdan et al. (2003) 
is conhned to a region about 400 km wide within the hux element. In our simulations, the 
transverse motion of the entire magnetic element creates strong pressure perturbations just 
outside the tube which in turn are responsible for the vortical motions and the generation 
of acoustic waves in the ambient medium as well as the strong longitudinal motions that 
eventually produce shock heating in the upper atmosphere. The recognition that the hux 
tube interface is a source of acoustic waves in the ambient medium and that these waves 
can penetrate back into the tube at higher elevations is a new feature of the present work. 
Thirdly, in addition to periodic driving, we also consider impulsive motions at the lower 
boundary. 

We should point out some limitations of our current study. Our analysis is based on a 
two-dimensional treatment in slab geometry, assuming that huid displacements are conhned 
to the X — z plane. The waves examined by us are diherent from the kink and sausage modes 
treated earlier by us using a simplihed thin hux tube approximation. Ideally, one would 
use a 3-D treatment which is presently beyond the scope of this work. Our analysis also 
neglects the torsional Alfven wave, which in reality would couple to the fast and slow modes. 
Even within the two-dimensional framework, we have had to restrict the height range in our 
simulations to 1200 km, since the Alfven speed (that essentially controls the time step used 
in the numerical scheme) increases very rapidly with ; 2 . 
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In summary this paper is the hrst of a series of investigations on the dynamics of the 
solar magnetic network based on a two-dimensional treatment of the MHD equations. We 
have made a start by considering processes occurring in a single flux tube which expands 
with height in the photosphere and assumes a ’’wine glass” geometry in the chromosphere. 
Flows and wave motions are generated in this conhguration by transverse motions at the 
base of the flux tube. For impulsive driving we find the presence of vortical motions. An 
interesting feature for both impulsive and periodic driving, is the development of shock-like 
features in the upper atmosphere which can be important in heating the chromosphere. We 
hope to extend the scope of the present calculation in future work to include multiple flux 
tubes as well as include a larger height extension in order to estimate transport of energy in 
to the corona through this mechanism. 


SSH thankfully acknowledges support from the National Science Foundation, U.S.A. 
through grant number ATM-0207641. 
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Table 1. 


Parameters on the Tube Axis & Ambient Medium 


Variable 


Tube Axis 

Ambient Medium 

Base 


Top 

Base 

Top 

Temperature 

4800 K 


8200 K 

4800 K 

8200 K 

Density 

1.3 10 ^ g cm ^ 

9.5 10“^^ g cm“^ 

4.0 10“^ g cm“^ 

2.9 10“^^ g cm“^ 

Pressure 

4.1 10^ dyn 

cm ^ 

8.3 dyn cm“^ 

1.2 10^ dyn cm“^ 

22 dyn cm“^ 

Sound speed 

7.1 km s“^ 


9.3 km s“^ 

7.1 km s“^ 

9.3 km s“^ 

Alfven speed 

11 km s“^ 


93 km s“^ 



Magnetic field 

1420 G 


102 G 



/9 

0.5 K 


0.02 
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RADIUS (Mm) 

Fig. 1.— Model of a network element consisting of individual flux tubes separated at the 
photospheric surface by a distance of 1000 km which merge at a height of about 600 km. 
The red box corresponds to the domain taken up for dynamical simulations. 
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Fig. 2.— Variation of (a) the temperature T, and (b) the logarithm of the pressure p as a 
function of height on the held lines corresponding to / = 0.55 (solid curve), 0.75 (dotted 
curve) and 0.95 (dashed curve) in the equilibrium conhguration. Note the temperature is 
uniform (by assumption) in the horizontal direction. 
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Fig. 3.— Variation of (a) the magnetic held strength B, and (b) /3 as a function of height 
on the held lines corresponding to / = 0.5 (solid curve), 0.8 (dotted curve) and 0.9 (dashed 
curve) in the equilibrium conhguration. 
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Fig. 4.— Variation of with horizontal distance (from the left edge of the computational 
domain) at the following heights: 2 ; =0 (black solid line), 500 km (black dotted line) and 
1000 km (black dashed line). The corresponding red curves denote measured with respect 
to the values given on the right scale. 
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Fig. 5.— Flow pattern and temperature perturbation, , AT, (about the initial state) at (a) 
55 s, (b) 82 s, (c) 109 s and (d) 136 s in a network element due to horizontal motion at the 
lower boundary, with an amplitude of 750 m s“^, applied for half a wave period {P = 240 s), 
after which the motion ceases. The thin black curves are the magnetic held lines, the arrows 
denote the direction of the how, and the color scale shows the temperature perturbation. 
The white curves denote contours of constant [3 corresponding to /? = O.l, 0.5, l.O (thick 
curve) and lO. 
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Fig. 6.— The relative pressure perturbation (about the initial state), Ap/po, at t = 82 s in 
a network element due to horizontal motion at the lower boundary with an amplitude of 750 
m s“^ and a period of P = 240 s. The thin black curves are the magnetic held lines and the 
white curves denote contours of constant f3 corresponding to {3 = 0.1, 0.5, 1.0 (thick curve) 
and 10. 
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Fig. 7.— Variation of the velocity components (a) along the magnetic held Vg and (b) 
normal to the magnetic held Vn at diherent t = 82 s (solid curve), 109 s (dotted curve) and 
150 s (dashed curve) on a held line characterized by / = 0.6 in a network element due to 
a horizontal motion at the lower boundary with an amplitude of 750 m s“^ and a period 
P = 240 s. 
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Fig. 8.— Flow pattern (arrows) and the temperature perturbation (about the initial state), 
AT, at (a) 41 s, (b) 68 s, (c) 96 s and (d) 123 s in a network element due to a periodic 
horizontal motion at the lower boundary with an amplitude of 750 m s“^ and a wave period 
of P = 24 s. The thin black curves are the magnetic held lines. The white curves denote 
contours of constant f3 corresponding to {3 = 0.1, 0.5, 1.0 (thick curve) and 10. 
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Fig. 9.— The normal velocity component at t = 68 s in a network element due to a 
horizontal motion at the lower boundary with an amplitude of 750 m and a period of 
F = 24 s. The thin black curves are the magnetic held lines and the white curves denote 
contours of constant /3, corresponding to (3 = 0.1, 0.5, 1.0 (thick curve) and 10. 
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Fig. 10.— Variation of the velocity component along the magnetic field, Vg, at times t = 68 s 
(solid curve), 96 s (dotted curve) and 123 s (dashed curve) on a held line characterized 
by / = 0.6 in a network element due to horizontal motion at the lower boundary with an 
amplitude of 750 m s“^ and a period of P = 240 s. 











